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We analyze the non-local transport properties of a d-wave superconductor coupled to metallic 
electrodes at nanoscale distances. We show that the non-local conductance exhibits an algebraical 
decay with distance rather than the exponential behavior which is found in conventional supercon- 
ductors. Crossed Andreev processes, associated with electronic entanglement, are favored for certain 
orientations of the symmetry axes of the superconductor with respect to the leads. These properties 
would allow its experimental detection using present technologies. 

PACS numbers: 74.20.Rp,74.50.+r,74.45.+c,81.07.Lk 



I. INTRODUCTION 



Cooper pairs in superconducting nanostructures pro- 
vide a potential source of entangled electrons [J, @, H, 
EL HI, a possibility that has been recently explored in 
conventional superconductors both theoretically [1, 0, 
i^JTSJlljJlJ, O 0, E [3 and experimentally 
[17, Il8l 1 19L 2CJ|. In a typical experimental device, a 
superconducting region is contacted by several metallic 
electrodes at nanoscale distances with the aim of analyz- 
ing the non-local transport properties at subgap voltages. 
In the limit of vanishing contact transparency the non- 
local conductance is controlled by two type of processes 
yielding opposite contributions: direct elastic tunneling 
of electrons between two separate leads (elastic cotunnel- 
ing, EC) and crossed Andreev reflection (CAR) processes 
in which injected electrons from one lead are reflected as 
holes in the other lead (see Fig. [I]). The time reverse of 
these last processes involve entangled electron pairs on 
two separate leads [2l[ . In conventional superconductors 
the average conductance tends to cancel due to the oppo- 
site contribution of EC and CAR processes . Several 
mechanisms have been proposed to avoid a complete can- 
cellation and have been invoked to explain the available 
experimental results. Among them one can quote the 
use of ferromagnetic leads [J EH EH E2 and the effect 
of electron-electron interactions in certain experimental 
geometries [lSj . On the other hand, the magnitude of 
these non-local processes decays exponentially with the 
distance between the leads on a scale fixed by the super- 
conducting coherence length £o- In practice this means 
that non-local effects can be observed on distances of the 
order 10 nm to 1/im, depending on the material (l7l. [Tc|. 

These effects have been much less explored in the case 
of unconventional high critical temperature superconduc- 
tors (HTcS) [6l.l22l.l23l]. The characteristic small values of 
the coherence length in these systems cast doubts about 
the observability of non-local correlations. However, due 
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FIG. 1: Two leads on a d-wave superconductor. The distance 
between the leads is d, a is the angle between the crystalo- 
graphic axes of the superconductor and the direction normal 
to the surface. An incoming electron from lead a can be re- 
flected in lead b as: (b) an electron or, (c) a hole while a 
Cooper pair is created in the superconductor. 



to the anisotropy of the pair potential, the coherence 
length along certain directions can be much larger than 
£o- This anisotropy is behind the non-local nature of the 
electromagnetic response of HTcS[24|, [25[. In fact, some 
indirect evidence of CAR processes in HTcS coupled to 
ferromagnetic leads has been presented [2(| [27| • 

In the present work we analyze the non-local transport 
in d-wave superconductors and show that in contrast to 
the conventional s-wave case, CAR processes are long 
ranged. Moreover, we show that for certain orientations 
of the axes of the superconductor with respect to the con- 
tacts, CAR processes dominate over EC at low voltages 
and small contact transparency. We believe that these 
findings open the possibility of using HTcS as a source 
of entangled electron pairs. 



II. CROSSED ANDREEV REFLECTION AND 
ELASTIC COTUNNELING IN HIGH T c 
SUPERCONDUCTORS 



The situation to be analyzed is illustrated in Fig. Q] 
We consider a semi-infinite d-wave superconducting re- 
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gion connected to two normal leads, denoted by a and 6, 
and separated a distance d. Our aim is finding the cur- 
rent induced on lead b, lb, when a voltage V a is applied 
on lead a. The two processes contributing to this current 
are depicted on panels (b) and (c). Being the result of the 
diffraction of quasiparticles by an anisotropic pair poten- 
tial, the relative weight of the two processes will be af- 
fected by the orientation of the superconductor symmetry 
axis (angle a in Fig. [1]). In fact, EC processes are favored 
for electron propagation along the nodal lines (where the 
order parameter vanishes), while CAR processes reach a 
maximum amplitude along directions where the modulus 
of the pair potential is maximum. Although on a spatial 
average the contribution of the two processes should be 
equal as in the case of isotropic s-wave superconductors, 
a dominance of one of the two can be found for spe- 
cific orientations of the symmetry axis with respect to 
the leads. These qualitative arguments allow to under- 
stand the dominance of CAR over EC processes for the 
a = case [d x 2__ y 2 symmetry) and the opposite behavior 
in the case of a = 7r/4 (d xy symmetry) as discussed in 
detail below. 

In the spirit of the Hamiltonian approach of Ref. [28[ 
the differential conductance defined as Oba = dIb/dV a , 
can be written as [l(| 



a ba = 8 " 2 ^ Pe,a(p M |GL, 12 (eF)| 2 -Pe,^ a|11 (eV)| 2 ), 

(1) 

where p e (h),a{b) is the local density of states of the elec- 
tron(hole) in the lead a(b), while p a and pb denote the 
corresponding hopping parameters coupling the super- 
conductor to the leads. The quantities G£ a n (eV) and 
G r ha 12 (eV) are the non-local propagators in the super- 
conducting region (indexes 1 , 2 refer to electrons and 
holes in Nambu space). The first term on the right 
hand size of Eq. (TTJ) is due to CAR processes, while the 
second term corresponds to EC. Therefore the crossed 
differential conductance is positive if CAR dominates 
over EC or negative in the opposite case. In order to 
make contact with possible experiments it also conve- 
nient to analyze the non-local resistance Rb a , given by 
-Vba/{(TaaVbb - Vab&ba), where er QQ ( b& ) are the local con- 
ductances which can be obtained within the same for- 
malism [T3, UK ■ The propagators G r a p ~ of the coupled 

system are then given by G r (E) = (g r (i?) _1 + it) , 
where g r is the retarded Green function of the uncoupled 



superconductor and T a fi,ij — PairpN8 a ,p§i,j' We have 
assumed that the densities of states of the normal metals 
are energy independent, i.e. p e , a (b) = Ph,a(b) = PN- The 
symbol v here denotes 4x4 matrices defined in the elec- 
trodes ffi Nambu space, while we reserve the symbol A for 
the reduced 2x2 Nambu space. To calculate g we first 
determine the superconductor surface Green function in 
momentum representation, gg{E,k y ), using the asymp- 
totic solutions of the Bogoliubov de Gennes equation [2!| , 
which yields 
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k± = ^k 2 xF ±2mn/h 2 ,T=\A+\/(E + n) 
<P± = arg(A ± ), Aip = ip + - tp_, k 2 xF = k 2 F -k 2 
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(1 - T 2 e- iA f) , n = ^E 2 - |A+| 2 
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In the above equations A is the pair potential which 
depends on the wave vector, taking the values A + and 
A_ along the directions and ir — 8 respectively, where 
8 = tan -1 {ky/kxp)- These are given by A±(0) = Ao, for 
s symmetry and A± = Ao cos(2 (8 T a)) for (i-symmetry. 
The retarded component is obtained by adding a small 
positive imaginary part irj to the energy. From gs(E, k y ) 
one then obtains the non-local components g r ha by 
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g r s {E,k y )\f(k y )\ 



dk 



(4) 



where the weighting factor f(k y ), proportional to the 
perpendicular wave vector k x p, provides the appropri- 
ate connection between the continuous model used to 
describe the superconducting region and the discrete 
Hamiltonian approach used to obtain Eq. |T]) (see Refs. 

dill). 

As a first test of the model one can check that in 
the case of s-symmetry for E < Aq and kpd >> 1, 



cos (k F d) - sin 2 (k F d)) /d 3 with £(E) 



aba « e 

£o/Re(v/l - E 2 /A 2 ) andfo = ^/(ttA ), a result 
which agrees with Refs. [OiEl- Therefore aba exhibits 
changes in sign on the A p scale and its spatial average is 
zero la. fl2l. 



A. Results for d x 2_ y 2 symmetry 

We now consider the d x 2_ y 2 symmetry. Due to the 
anisotropy of the pair potential an incoming electron 
from lead a is scattered as a quasiparticle in the su- 
perconductor, exploring regions where E > A(0) and 
E < A(0), with an effective coherence length £(E,8) = 
£o/Re(yI — E 2 /A(8) 2 ) which takes values from £o to oo. 
For this reason one typically finds that the propagators 
exhibit a slower decay with distance than in the case of 
s-symmctry. In this paper we have fixed Ao ~ 20meV 
Ml and A /Ep ~ KT 1 as typical values for HTcS [13]. 
We also take rj ~ 0.002A to simulate the effect of weak 
disorder [33[. Figure [2] illustrates the spatial dependence 
of the Green functions. Due to the dependence on k y of 





FIG. 2: (Color online) The anomalous propagator |<?i, a ,i2| 2 for 
d x 2_ y 2 symmetry at zero energy as a function of the separa- 
tion between the leads. The propagator has been normalized 
to its value in the normal state at d = 0. The correspond- 
ing electron propagator gb a ,n is negligible within this scale. 
The envelope curve, decaying as 1/d 2 is indicated by a full 
line. The corresponding curve for s-symmetry with the same 
choice of parameters, exhibiting an exponential decay, is also 
represented for comparison. The inset shows the exponents 
n,2 in the decay laws of Eq. ([5]) as a function of energy. 

the pair potential, it is not possible to obtain an ana- 
lytical expression of their variation with d as in the case 
of s symmetry. However, from numerical regressions for 
kpd » 1, |.9f, a ,n| 2 and g& a ,i2| 2 can be fitted as 

\„r (ww 2 ~ c m+ d m cos 2 (fcd) 

|S6a,H(12)WI " \ kFd \r im ' W 

The values of the exponents f 1(2) fixing the spatial decay 
are shown in the inset of Figj2] as a function of energy. 
For low energies (E « A ) k ~ kp, cui) « dip) and 
di » d\ — ► for E — * 0, and therefore the propagator 
\dba 12 1 takes a much larger value than \g£ a n |, yielding 
a clear dominance of CAR over EC in the tunnel limit. 
Notice that the low energy excitations at 9 ~ 7r/4 give 
a negligible contribution to \g£ a n | due to the weighting 
factor in Eq. (j4|) that is maximum at low angles. In 
contrast, most of the weight in \gl a 12 | comes from 9^0 
where A reaches a maximum. On the other hand, for 
energies higher than E ~ 0.1 A \gl a 12 | and \gl a 11 | tend 
to have the same magnitude on average. 

Fig. [3]further illustrates the different behavior of CAR 
and EC contributions to the non-local conductance for 
this orientation as one of the leads moves inside the su- 
perconductor while the second is located at the surface. 
These maps clearly correspond to a difraction pattern for 
electrons injected at one point in the surface. In spite of 
its complex structure one can identify the region of low 



FIG. 3: (Color online) Plots of the CAR and EC contribu- 
tions to the non-local conductance for the d x i_ y i symmetry 
as one of the contacts moves inside the superconductor while 
the other remains fixed at x — y = 0. The contacts are in 
the tunnel limit and the voltage is set to zero. Both contri- 
butions are plotted in a logarithmic scale normalized to their 
maximum value. 




FIG. 4: (Color online) Spatial averaged non-local resistance 
at d = 10£o for d x 2_ y 2 symmetry as a function of the voltage 
for different values of the transmission Tjv. The inset shows 
the non-local conductance at eV — as function of d and for 
different values of Tjv. For Tjv ~ 0.9 there is a sign change in 

Cba- 



angles from the surface (tt/2 > 9 > 7r/4) where CAR 
processes have a clear dominance and the nodal lines 
(9 ~ 7r/4) around which EC processes are favoured. 

The results for a^a and Rb a in the d x 2_ y 2 orientation 
for arbitrary contact transmission are illustrated in Fig. 
|4] It is found that an increase in transmission leads to 
a reduction of the CAR contribution while the EC one 
increases. As a consequence both quantities exhibit a 
change of sign when the coupling to the leads increases. 
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FIG. 5: (Color online) Non-local conductance at eV = for 
dxy symmetry as a function of the distance between the leads 
for Tat = 0.04,0.15 and 0.3. The inset shows the spatial 
averaged non-local conductance at d — 10£o and Tjv = 0.04 
as function of the voltage. 



This is illustrated for aba in inset of Fig. |4j We can 
obtain further insight on this effect at low energies where 
Eq. ([T]) can be approximated as 
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(6) 

In obtaining this expression we have assumed symmet- 
rical contacts (p a = Pb = p) with P = pup^ being 
the normalized hopping parameter, such that the normal 
transmission for a single contact is Tjy = 4P 2 /(1 + P 2 ) 2 . 
Within this approximation the dependence with the sepa- 
ration between the leads does not change when the trans- 
mission is increased, as it is seen in the inset of Fig. 2] 
On the other hand, this equation predicts a change in 
sign of aba for P ~ 0.72 (T/v — 0.9) in agreement with 
the numerical results in the inset of FigHJ The non-local 
resistance Rb a averaged on a range ~ Xf is shown in Fig. 
[4] for d = 10£o- We observe that in the low transmission 
regime this quantity is negative, as it corresponds to the 
dominance of CAR processes, and exhibits a peak at low 
bias. 



with al/E dependence in g£ a n . On the other hand, the 
distance dependence for low energy and kpd >> 1 is in 
this case approximately 1 /d 4 both for CAR and EC pro- 
cesses. In this orientation, the local Andreev reflection is 
zero because of diffraction of quasiparticles in the contact 
[35L f36j| and the CAR contribution to ab a is not zero, but 
is always smaller than the EC one, leading to a negative 
non-local conductance as shown in Fig. [5] Basically, the 
dominance of the EC contribution is caused by the sup- 
pression of the pair potential along the 8 = line. The 
effect of varying the contact transmission can be under- 
stood analytically within a similar approximation as done 
for the d x 2_ y 2 case, which allows to obtain the following 
expression for the crossed differential conductance 
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(7) 



Notice that CAR and EC contributions are equally af- 
fected by the coupling to the leads (through the P- 
dependent common denominator) and therefore EC dom- 
inates over CAR for the whole transmission range. The 
spatially averaged non-local conductance is negative and 
presents a zero bias peak that decreases with increasing 
transmission. 



III. CONCLUSIONS 

In summary we have analyzed the behavior of the 
crossed differential conductance in <i-wave superconduc- 
tors in a multiterminal configuration. We have shown 
that correlations between different leads exhibit an alge- 
braical decay instead of the exponential behavior which is 
typically found in conventional superconductors. In the 
case of d x 2_ y 2 orientation crossed Andreev processes are 
favored at low voltages and contact transmissions, while 
for the d xy case a zero bias non-local conductance peak 
appears, which is dominated by elastic-cotunncling. In 
both cases the spatially averaged non-local conductance 
is different from zero. These properties would allow to de- 
tect non-local transport at distances several times larger 
than the characteristic coherence length in these systems. 
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B. Results for d xy symmetry 

The results for d^-symmetry (a = 7r/4) are shown in 
Fig. [5l The main effect for this s ymm etry is the appear- 
ance of a zero energy bound state [34[ , which is associated 
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